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1. Introduction 

Whenever results obtained from lattice QCD simulations are to be confronted with experi- 
mental results it is important to have a sound control of systematic uncertainties emerging in lattice 
QCD. The most prominent of those are discretisation errors, finite size effects (FSE) and uncertain- 
ties arising from the unphysically large mass values usually simulated. The main reason for their 
prominence is the fact that lattice QCD simulations become increasingly computer time demanding 
when a) the lattice spacing is reduced b) the quark masses are reduced towards the physical point 
and c) the volume is increased. 

The control of these systematic uncertainties requires simulations with an ff{a) improved 
lattice formulation at sufficiently small values of the lattice spacing a, say a < 0.1 fm where &{a^) 
lattice artifacts are small. Physical volumes should be large enough, say with spatial box size L 
larger than 2 fm and mps •L>2> (wips is the mass of the lightest pseudo scalar particle). And, 
in order to be able to utihse chiral perturbation theory (;tPT) to bridge between simulated quark 
masses and the physical point, simulations with a range of masses are needed, with the smallest 
value of /Mps < 300 MeV. It goes without saying that the aforementioned bounds are only estimates 
and need to be checked carefully in actual simulations. 

Due to recent algorithmic improvements ijl]-^ it became possible to meet all these require- 
ments using Wilson's original formulation of lattice QCD. It has the advantage of being conceptu- 
ally clear and simple. And, i?'{a) improvement can be implemented in several ways, one of which 
is using so called Wilson twisted mass fermions at maximal twist. As was shown in Ref.[||], 
in maximally twisted mass lattice QCD (Mtm-LQCD) physical observables can be obtained ^(a) 
improved by tuning a single parameter only. In particular, no operator specific improvement co- 
efficients need to be computed. This theoretical expectation could be verified in the quenched 



approximation to work very well [g|-|12|] (for a recent review see Ref. []13||.) 

Based on these successes in the quenched approximation the European Twisted Mass (ETM) 
collaboration decided to start a large scale simulation project using two flavours of mass degenerate 
quarks with the lattice formulation of Mtm-LQCD. First accounts of this effort are published in 



Refs. [|14]-|16|] indicating that 0'{a) improvement works very well when sea quark effects are taken 
into account in the simulations. This proceeding contribution aims to summarise the progress and 
current status of the two flavour project of the ETM collaboration. 

2. Gauge and Fermionic Action 

In the gauge sector we employ the so-called tree-level Symanzik improved gauge action (tl- 



Sym) llnp, viz. 



bo X {l-^^Tv{Ull\)} + b, X {l-^^Tv{Ulll)}\ , 



\ 



with the bare inverse gauge coupling j8, bi = —1/12 and Z?o = 1 — SZji. The fermionic action for 
two flavours of maximally twisted, mass degenerate quarks in the so called twisted basis [^, ^ 
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reads 

5tm = [D[U]+mo + WT^] zW} , (2.1) 

X 

where mo is the untwisted bare quark mass, pLq is the bare twisted quark mass, is the third PauU 
matrix acting in flavour space and 

D[t/] = ^[7M(V;, + V;)-aV*V,] 

is the mass-less Wilson-Dirac operator. and V* are the forward and backward covariant dif- 
ference operators, respectively. Twisted mass fermions are said to be at maximal twist if the bare 
untwisted mass mo is tuned to its critical value merit, the situation we shall be interested in. For con- 
venience we define the hopping parameter K= 1/(8 + 2amo). Note that we shall use the twisted 
basis throughout this paper. 

Maximally twisted mass fermions share most of their properties with Wilson's originally pro- 
posed formulation, but provide important advantages: the spectrum of 75 {D[U]+mQ + ipCqJ^) ■ (h.c.) 
is bounded from below, which was the original reason to consider twisted mass fermions. The 
twisted mass jXq is related directly to the physical quark mass and renormalises multiplicatively 
only. Many mixings under renormalisation are greatly simplified. And - most importantly - as was 
first shown in Ref. ^ physical observables are automatically ^(a) improved without the need to 
determine any operator-specific improvement coefficients. 

The main drawback of maximally twisted mass fermions is that flavour symmetry is broken 
explicitly at finite value of the lattice spacing, which amounts to ff{a^) effects in physical observ- 
ables. However, it turns out that this is presumably only relevant for the mass of the neutral pseudo 
scalar meson (and closely related quantities). We shall discuss this issue in more detail later on. 

2.1 0'{a) Improvement in Practice 

It is well established in the literature that ^ {a) improvement works very well in practice in the 
quenched approximation [^-[Tl|, |T8[ [l9|, [T^]. As an example we show in figure |l(a)| the essentially 
flat continuum exti^apolation in a^ of the pseudo scalar decay constant /ps in physical units (ro = 
0.5 fm was used to set the scale) for three reference values of the pseudo scalar mass as obtained 



in Ref. |11]. 



So far we did not discuss how maximal twist can be achieved in practise and there are various 



solutions to this problem. Emerging from the proof of 0'{a) improvement at maximal twist pOl , |2l[ ] 
the general prescription is to choose a parity odd operator O and determine amcdt such that O 
has vanishing expectation value at fixed physical situation for all lattice spacings. The physical 
situation can be fixed for instance by keeping mps in physical units fixed, where mps is the mass of 
the lightest charged pseudo scalar particle. One possible quantity to tune is the PCAC quark mass, 
defined as 

z.{doA",{x,t)p-m 

mpcAC= 2i:,(P«(x)P«(0)) ' ^ = ^'^' ^^-^^ 

where and P" are the axial vector current and the pseudo scalar density, respectively. 
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Figure 1: (a) /ps [MeV] as a function of fl^[fm^] in the quenched approximation [|rT]]. (b) Monte Carlo 
history oiP— (P) and AH for ensemble Ci . 



Tuning the PCAC mass to zero has been successful in the context of the aforementioned quenched 
investigations, in agreement with theoretical considerations [22, 23, 2C]. The collaboration follows 
the strategy to determine the value of amcdt at the lowest available value of ajJ. <C aAqcd. which is 
sufficient to guarantee iff (a) improvement [pO[]. 



3. Setup and Tuning 

In table |l] we summarise the various ensembles produced by the ETM collaboration. We 
have simulations for three different values of the inverse gauge coupling j8 = 3.8, jS = 3.9 and 
j8 = 4.05. The corresponding values of the lattice spacing are about a w 0.10 fm, a f» 0.09 fm and 
a ^ 0.07 fm, respectively (see later). For each value of j8 we have four or five different values 
of the bare twisted mass parameter ajJ^q, chosen such that the simulations cover a range of pseudo 
scalar masses between 300 and 700 MeV. 

The physical box sizes of the simulations at j3 = 3.9 and j3 = 4.05 are roughly equal and 
around 2.2 fm, while the volume at j8 = 3.8 is slightly larger. For all j8 -values we have carried out 
simulations at different physical volumes in order to check for finite size (FS) effects. For each 
value of j3 and /i^ (ensemble) we have produced around 5000 equilibrated trajectories in units of 
T = 0.5. The actual trajectory length z is given in table In all cases we allowed for at least 1500 
trajectories for equilibration (again in units of T = 0.5). 

Note that the analyses for the ensembles at j3 = 3.8 are in a very preliminary status for reasons 
that will be explained later. For the purpose of this proceeding contribution j3 = 3.8 results are 
excluded from most of the present analyses. 

The determination of for instance quark masses requires a renormalisation procedure. To 
this end we have implemented the non-perturbative RI-MOM renormalisation scheme [p4|], which 
provides in the case of maximally twisted mass fermions ff{a) improved renormalisation constants 
[p5|]. Where appropriate, we convert our results to the MS scheme at the desired scale using 
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Ensemble 


L^xT 


P 




K 


Tint(P) 


Tint (amps) 


T 


Ai 


24^ X 48 


3.8 


0.0060 


0.164111 


190(44) 


8(2) 


1.0 


A2 






0.0080 




172(80) 


10(2) 


1.0 


A3 






0.0110 




130(50) 


6(1) 


1.0 


A4 






0.0165 




40(12) 


6(1) 


1.0 


A5 


20^ X 48 


3.8 


0.0060 


0.164111 


250(100) 


5(1) 


1.0 


Si 


243 X 48 


3.9 


0.0040 


0.160856 


47(15) 


7(1) 


0.5 


Bi 






0.0064 




23(7) 


17(4) 


0.5 


B3 






0.0085 




13(3) 


10(2) 


0.5 


B4 






0.0100 




15(4) 


7(2) 


0.5 


B5 






0.0150 




30(8) 


20(6) 


0.5 


Be 


32^ X 64 


3.9 


0.0040 


0.160856 


37(11) 


2.8(3) 


0.5 


Ci 


32^ X 64 


4.05 


0.003 


0.157010 


18(4) 


7(1) 


0.5 


C2 






0.006 




10(2) 


9(2) 


0.5 


C3 






0.008 




13(3) 


7(1) 


0.5 


C4 






0.012 




5(1) 


4.8(6) 


0.5 


Cs 


243 X 48 


4.05 


0.006 


0.157010 


12(2) 


11(1) 


1.0 


Ce 


20^ X 48 


4.05 


0.006 


0.157010 


10(2) 


7(1) 


1.0 



Table 1: Summary of ensembles produced by the ETM collaboration. We give the lattice volume x T 
and the values of the inverse coupling )3, the twisted mass parameter aHy, the hopping parameter K and 
the trajectory length T. In addition we provide values for the integrated autocorrelation time of two typical 
quantities, the plaquette P and the pseudo scalar mass amps, in units of T = 0.5. 



renormalisation group improved continuum perturbation theory at N LO [|26|]. For details on this 



procedure we refer to Ref. [ |25| , [27| |. 



3.1 Algorithm Stability and Performance 



The algorithm we used is a HMC algorithm [28] with mass preconditioning [|l], ^ and mul- 
tiple time scale integration (mtmHMC), as described in detail in Refs. [§, 30|. The algorithm per- 
forms smoothly and without any instabilities in the whole range of simulation parameters we have 
available. In particular we do not observe any instabilities at our smallest values for the twisted 



mass parameter ajx^ for the different values of j8. In figure |l(b)| we illustrate this with the Monte 
Carlo history of the plaquette P — (P) and the difference in the HMC Hamiltonian AH for ensemble 
Ci- For further details about the algorithmic parameters used see Ref. 

Due to the fact that the twisted mass parameter provides an infra-red cut-off for the eigenvalue 
spectrum of the lattice Dirac operator we do not expect instabilities due to very small or zero 



eigenvalues (for a discussion of this issue for Wilson fermions see Ref. [|32|]). However, Wilson 
type fermions exhibit a non-trivial phase structure at finite value of the lattice spacing with a first 
order phase transition at the chiral point for our choice of the gauge and fermionic action [^-^]. 
As a consequence simulations at given, finite value of the lattice spacing cannot be performed with 
arbitrarily small values of the bare quark mass and one needs to make sure that the values of the 
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Figure 2: (a) Comparison of the cost estimate we measure for the mtmHMC (data points) to the DD-HMC 
(hnes). For the DD-HMC we used Eq.(3.1). The upper (lower) hne compares to the j3 = 4.05 (j3 = 3.9) data 
points, (b) ro/a as a function of (ajJ.)^ for j3 = 3.9. The hne represents a linear extrapolation in {afJ.)^ to the 
chiral Umit. 



twisted mass parameter are large enough for the simulations not to be affected by the first order 
phase transition. 

This issue was thoroughly investigated by the ETM collaboration with the result that we do 
not see metastabilities for any of our simulation points. But, as a matter of fact, simulations at large 
value of the lattice spacing and maximal twist are performed potentially in the close vicinity of the 
(second order) critical endpoint at ajJ-c of the first order phase transition line. This line extends in 
the /x^-K"-parameter plane to a first approximation perpendicular to the K"-axis from (K^rit) — Mc) to 
{Kcrit, IJ-c) (see Ref. [ |T3| ] and references therein for a more detailed discussion). Hence, there is the 
danger for long autocoiTclation times in quantities that are not continuous at the phase transition 
like for instance the plaquette or the PC AC quark mass mpcAC- 

This is supported by table [l| where we give estimates of the integrated autocorrelation times for 
the plaquette and the pseudo scalar mass. While for the pseudo scalar mass Tim is always moderately 
small and depending only weakly on the values of j8 and /x^, this is not the case for the plaquette: 
there is a trend in the data that Tint{P) decreases with increasing values of /i^ and j8, even though 
the statistical errors are so large that the /^^^ -dependence is not significant. Note that the Tint-values 
for nipcAC are similar to those for the plaquette. Fortunately, we observe very long autocorrelation 
times only for our smallest j8 -value. This is the reason that we are stiU investigating the error 
analysis for those ensembles. This affects the determination of the PCAC quark mass, which is 
needed to tune to maximal twist. Hence, at this stage, we use the corresponding results at j8 =3.8 
only for estimating systematic errors. 

In order to investigate algorithm stability it was suggested in Ref. [ p9| ] to study the statistical 
distribution of the twist angle as a function of the quark mass. The quantity of interest is the PCAC 
quark mass, evaluated in the background of a given gauge configuration by taking the axial to 
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pseudo scalar correlator at some large value for f/a normalised by fl'mps/(2Cpp). Here Cpp is the 
pseudo scalar correlator averaged over all gauges evaluated at the same value oft/a. The PC AC 
quark mass is related to the twist angle ft) by 

tan ftj — ^ 



Zawipcac 



The statistical distribution of ampcAC resembles always to a very good approximation a normal 
distribution. As will be discussed later, at fixed value of j3and L/a the distribution mean depends 
on aiJ.q as expected. For all three j3 -values the standard deviation is increasing at the order of a few 
percent with increasing a/i^. In addition we observe a smooth dependence of the standard deviation 
on L/a and the lattice spacing. Hence, there is no sign for instabilities seen in tanftj. 

It is interesting to compare the performance of our HMC variant with the one using domain de- 
composition as a preconditioner (DD-HMC) (instead of mass preconditioning), which is described 



in Ref. [|40|]. A useful performance figure is the number of floating point operations Cop required 
to generate 1000 independent gauge configurations as a function of the box size L[fm], the lattice 
spacing ^[fm] and the quark mass m[MeV] in the MS scheme at = 2 GeV [^, ^ 

/20MeV\'^'" / L fOAfmY" 



In Ref. p2| ] the parameters in Eq. (3.1) for the DD-HMC algorithm with Wilson fermions are 
estimated roughly to k = 0.3, c„j = 1, cl = 5 and Ca = 6, which is a significant improvement as 
compared to cost estimates for the original HMC algorithm, see for instance Ref. pi]]. Using 
the integrated autocorrelation time of the plaquette as a measure for the autocorrelation of two 
gauge configurations we can measure Cop for the mtmHMC algorithm and compare the result in 



figure 2(a)| to the cost estimate for the DD-HMC algorithm [ i2 1 



Figure 2(a)| reveals that the performance of the two algorithms is very similar. In particular 



for the larger value of the two plotted lattice spacings the agreement is rather good. For our HMC 
version the lattice spacing dependence appears to be milder. However, one should keep in mind the 
large errors associated to this cost figure. Moreover, the result may depend significantly on how 
much effort is invested into tuning of algorithmic parameters. 

3.2 Sommer Parameter 

In order to be able to compare results at different values of the lattice spacing it is convenient 
to measure the hadronic scale ro/a [p3|]. It is defined via the force between static quarks at inter- 
mediate distance and can be measured to high accuracy in lattice QCD simulations. For details on 



how we measure ro/a we refer to Ref. []3l|]. 



In figure 2(b) we show ro/a as a function of {apLqY for j8 = 3.9. The mass dependence appears 
to be rather weak and a quadratic fit in ajiq to our data describes the data rather well. The results 
we obtain for the three jS-values in the chiral limit are ro/a = 4.46(3) at j8 = 3.8, r^/a = 5.22(2) at 
j8 = 3.9 and r^/a = 6.61(3) at j8 = 4.05. The statistical accuracy for ro/a is about 0.5%. Note that 
the data is also compatible with a linear dependence on apLq and a linear fit gives consistent results 
in the chiral limit. 
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Figure 3: (a) renormalised PCAC quark mass against renormalised twisted mass for /3 = 3.9 and j3 = 4.05. 
The statistical uncertainty on Zp is not included, (b) Effective mass plot for the charged pseudo scalar mass 
for ensemble By, B3 and B5. The lines represent the fitted value for amps. 



3.3 Tuning to Maximal Twist 

In order to obtain ^ {a) improvement the bare quark mass /mq must be tuned to its critical value. 
As mentioned before we use the PCAC quark mass defined in Eq. (22) for tuning to maximal 
twist. The goal was to tune this quantity to zero approximately at the smallest available /i^-value 
at each lattice spacing, which corresponds to approximately fixed physical pseudo scalar mass. 
Considering only j3 = 3.9 and j8 = 4.05, it was possible to perform this tuning task with two or 
three tuning runs for each lattice spacing. 

The result is plotted in figure [3(a)| , where we plot the renormalised mpcAC against the renor- 
malised pLq, both in units of tq, for j8 = 3.9 and j3 = 4.05. The renormalisation factors were 
determined using the RI-MOM scheme. Within the statistical accuracy the PCAC quark mass is 
zero at a common value of the renormalised twisted mass. For all other values of pLq we observe 
(small) deviations from zero. This pLq dependence is an ^{a) cut-off effect which will modify only 
the 0'{a^) lattice artifacts of physical observables. The numerical precision we were aiming for 
was cjAqcd^/M i$ 0.01, where e is the uncertainty on the PCAC quark mass at the lowest value of 
/i^ [ 44 ] . It turns out that the accuracy we have achieved is sufficient for excellent scaling properties, 
as will be seen later. 



4. Results 

The ETM collaboration is currently analysing the generated gauge configurations for many 
different observables. Not all of these - partly preliminary - results can be summarised here and 
we shall therefore mainly concentrate on the pseudo scalar sector. At the end of this section we 
present an overview of the available results and give references of the corresponding proceedings 
contributions. 
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Figure 4: (a) ro/ps as a function of (romps)^ for j3 ~ 3.9 (Bi to B5) and j3 = 4.05 (Ci to C4). (b) Continuum 
extrapolation of /ps at fixed volume for three reference values of rQmps- The data points at j3 = 3.8 are not 
used. 



4.1 Charged Pseudo Scalar Mass and Decay Constant 

The charged pseudo scalar meson mass amps can be extracted from the time exponential decay 



of suitable correlation functions. For details on our analysis procedure see Ref. [31 1. To demon 



strate the quality of our data we show in figure |3(b)| effective mass plots for the ensembles Bi, S3 



and B5 obtained from the pseudo scalar correlation function. The final values for the masses are 



obtained from a fit to a 4 x 4 matrix of correlators [31 1. We also attempted to determine the energy 
of the first excited state of the pseudo scalar meson. We were not able to determine it with any 
reliability from an unconstrained fit. (In particular, there is no evidence for an excited state with 



mass + mp^, which is theoretically possible for maximally twisted mass fermions [13].) Con- 
straining the energy of the first excited state to three times the ground state mass, however, does 
allow for an acceptable fit. 

For maximally twisted mass fermions the charged pseudo scalar decay constant a/ps can be 
extracted from 

M = ^mn^)\, a = \,2 (4.1) 

due to the exact lattice PCVC relation with no need to compute any renormalisation constant [^]. 
Thanks to this advantage and having the results for ro /a at hand we plot the results for ro/ps as 
a function of (romps )^ for j8 = 3.9 and j8 = 4.05 in figure 4(a) . We plot only the results for the 



ensembles Si to S5 and Ci to C4 in order to have approximately equal physical volumes for the 
two values of jS . Figure 4(a)| provides first evidence that lattice artifacts in those two quantities are 



small. 



Continuum Extrapolation of fps at Fixed Volume 



This statement can be brought to a more quantitative level for instance for the pseudo scalar 
decay constant. To this end we interpolate a/ps at every lattice spacing linearly in (awips)^ to fixed 
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Figure 5: (a) Relative finite size effects for ensembles C2, C5 and Cg. The line represents a fit with for- 
mula (4.3) to our data, (b) /ps//o as a function of X^i/fo obtained from a combined fit to ;|fPT formulae. 



The dashed line is the fitted function Eq. (p~7|), and the data points are FS corrected. 



values of romps. We corrected for the very small difference in the physical volume [Q]. 

The interpolated data points are plotted in figure ^(b) for three values of rowips = 0.70, r^mps = 
0.90 and rowips = 1.10 as functions of (a/ro)^. It is visible that the differences between the data at 
j8 = 3.9 and jS = 4.05 are of the order of the statistical accuracy and hence, we perform a weighted 
average of these two lattice spacings to obtain continuum estimates. This indicates again that 
scaling violations are very small in the pseudo scalar decay constant. 



Even though we are still investigating the error analysis we plot in figure 4(b) also the results 
for jS =3.8 for those reference points where we are able to interpolate to the reference values of 
romps (for the lowest value of romps we currently need to extrapolate). Keeping in mind that these 
data points are very preliminary we can nevertheless say that they fit rather well into the picture 
that lattice artifacts are very small. 

Other quantities - for instance the renormalised quark mass - show a similar scaling behaviour 
as discussed in Ref. [Q], where also the results at j3 = 3.8 are discussed in more detail. 

Finite Size Effects 

At the level of statistical accuracy we have achieved now, finite size effects for /ps and mps 
cannot be neglected. It is therefore of importance to study whether the FS effects can be described 
within the framework of chiral perturbation theory. This requires to compare simulations with 
different lattice volumes and all other parameters kept fixed, like for instance ensembles C2, C5 
and Cg or Bi and Bg. For all these ensembles mpsL > 3 holds, which is believed to be needed for 
XVT formulae to apply. The first observation we make is that the finite size effects are compatible 



with an exponential behaviour in mpsL. As an example we plot in figure p(a)| the relative finite size 
effects 

0(L) - Oil = 00) 
Ro = — ^- — - (4.2) 
^ 0{L = 00) 
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Ensemble 


mpsL 


^meas 




^CDH 


mps 


By 


3.3 


+1.8 


+0.62 


+ 1.0 


/ps 


Bi 


3.3 


-2.5 


-2.5 


-2.4 


nips 


Cs 


3.5 


+1.1 


+0.8 


+1.3 


/ps 


C5 


3.5 


-1.8 


-3.2 


-2.9 


mps 


Ce 


3.0 


+6.2 


+ 1.8 


+4.7 


/ps 


Ce 


3.0 


-10.7 


-7.3 


-8.9 



Table 2: Comparison of measured relative FS effects Eq. (4.2) in % to estimates from xPT^ formulae. 



for O = mps against mpsL. The value of mps = mps(L 
perturbation theory inspired formula 



00) was obtained by fitting the chiral 



mps(L) = mps + ae-'^Ps^/L^/^ 



(4.3) 



to our data, with a and mps being free parameters [)45|]. A similar fit can be performed for /ps by 
using the value of mps as obtained from the first fit as input. Both fits describe the data rather well 
with x^/<iof < 1. 

We shall now compare the measured finite size effects to predictions of continuum chiral per- 
turbation theory (;tPT) (lattice artifacts appear to be negligible, as discussed before). The compar- 
ison reveals that continuum ;^PT formulae can describe FS effects within our statistical accuracy. 

The NLO ;tPT formulae for nips and /ps were derived in Ref. [ p6| ] (for short GL) and can be 
written as 



1 



mps(L) = mps 1 + -^gi (A) = mps (L) , 



/ps(L)=/ps l-2^fi(A) =MKf^{L) 



(4.4) 



(4.5) 



where 

^ = mis/ {471 fof , A=mpsL, 

gi is a known function p6| ] and the finite size corrections ^T^j depend apart from L and mps 
only on the unknown leading order low energy constant /o (note that our normalisation is such 
that fjc = 130.7 MeV). For the pseudo scalar meson mass the corrections are also known to two 



loops [47], but the asymptotic Liischer formula [48, 49] (for short CDH) provides an easier way to 
access higher order corrections to mps and /ps, and the differences to the NNLO formula turn out to 



be small [147|]. The drawback of CDH compared to GL is that additional parameters are needed as 
an input, among others the low energy constants Ai, A2, A3 and A4. Since we do not have enough 
data points to determine all these parameters from a fit to only FS data we have to rely on estimates 



available in the literature [49] 



Assuming that the results for the ensembles Be and C2 provide the infinite volume estimates 
for mps and /ps, we can compute the relative FS effects using Eq. ( |4.2[ ) for the ensembles Bi , C5 and 
Cg, which we denote with 7?"^'=^^ These measured estimates can then be compared to the ;^PT pre- 
dictions computed with the formulae from GL and CDH, denoted by R'^^ and R'"^^, respectively. 
The values of the unknown low energy constants are set to the estimates provided in Ref. ||49|]. In 
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Figure 6: (a) (mps//o)^ as a function of Xnl fo (b) '"ps/Xm a function of x/fo- The data point for 



ensemble C5 is slightly displaced. In (a) and (b) the dashed line represents the fitted functions Eqs. (i.6) and 
(|4.7|), respectively, and the data points are FS corrected. 



order to do so and for evaluating the CDH formulae we need the value of the lattice spacing which 
we estimate using ro = 0.45 fm. 

The results for R'^^'^^^ R^^ and R^^^ are compiled in table ^ It turns out that - in particular for 
mps - the asymptotic formula from CDH describes the data better than the one loop formula from 
GL: the CDH corrected data is in agreement with the infinite volume estimate within the statistical 
accuracy. Note that the observation that GL usually strongly underestimates the FS effects in mps 



can also be made from the Wilson data published in Ref. [ ]50| ] (cf. Ref. []39|]). 

These results make us confident that our simulations have eventually reached a regime of 
pseudo scalar masses and lattice volumes where ;^PT formulae can be used to estimate FS effects. 
But it is clear that in particular the CDH formula is affected by large uncertainties, mainly stemming 
from the only poorly known low energy constants, which are needed as input. Changing their values 



in the range suggested in Ref. [i9\, however, changes the estimated finite size effects maximally at 



the order of about 20% (of the corrections themselves). 

Quark Mass Dependence and Chiral Perturbation Theory 

So far we have argued that lattice artifacts in charged nips and fps for j8 = 3.9 and j8 = 4.05 
are very small and that FS effects can be described using ;f PT formulae. We shall now show that 
also the quark mass dependence of these quantities can be successfully described using ;tPT. This 
will in addition allow the determination of most of the aforementioned low energy constants and 
the lattice scale. 

As chiral symmetry is broken by the lattice Wilson term, eventually the chiral extrapolation 
should be done using continuum extrapolated, infinite volume data. This is described for our data 



in Ref. [44|. However, since lattice artifacts are not visible with our current statistical accuracy of 



about 1 % (and chiral symmetry breaking as well as flavour symmetry breaking are formally lattice 
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artifacts of ff{a^)), we can follow a different approach. 

This approach consists of describing our data for /ps and mps at j3 = 3.9 and jS = 4.05 simul- 
taneously with continuum chiral perturbation theory. We fit the appropriate {Nf = 2) continuum 
NLO xVT formulae [|6[ ^ 



ml^{L,^,) = Kl{L) [1 + ^ \og{Xi,/K\)] , (4.6) 

/ps(L,At,) =/o Kf{L) [l-2^1ogU^/A2)] , (4.7) 

to our raw data for mps and fps simultaneously for j8 = 3.9 and j8 = 4.05. K,„j{L) parametrise the 
FS corrections for which we can either use the GL or the CDH formulae. Both and Kf depend 
on low energy constants as well as on L and mps. The notation is 

X^, = IBqZ^IIci , ^ = Xn/i^T^h? , (4.8) 
and the normahsation /o = ^/IFq, i.e. = 130.7 MeV. In Eqs. (O) and (O) NNLO ;tPT 



corrections are assumed to be negligible. This approach has the advantage that we can include 
finite size data consistently in the fit and that we can use more raw data points. In addition we 
do not need to interpolate our data to reference points. The fit presented in the following is an 
extension to the fit presented for only j3 = 3.9 in Ref. [|T^. 

The fit can be parametrised by six free parameters: two dimensionless ratios A3,4//o and /o 
and Z^Bo in lattice units for both j3-values, i.e. A3//0, A4//0, a/o|j3=3.9, '^/o|/3=4.05» '^^^i^o 1/3=3.9 
and aZ^ Bo 1/3 =4.05- 

Finite size effects are corrected for by using the asymptotic formulae from CDH, which is 
consistently included in the fit. The parameters that are not fitted, basically Ai and A2, are set 



to the values suggested in Ref. [)49[]. Also the lattice spacing in fm, needed for evaluating the 
CDH formula, can be determined consistently from the fit by setting /ps = /;;: where the ratio 
mps//ps assumes its physical value. Note that in this ratio we use the physical value of the neutral 
pseudo scalar meson mass, in order to account for electro magnetic effects not present in the lattice 
simulation. 

The ensembles we use in the fit are Bi to B4 and C\ to C3. In addition we include the ensembles 
and C5 in the fit in order to explore the L dependence of our data. We do not use Cg in order not 
to give too much weight to this /x^-value. The ensembles B=, and C4, the largest available masses, 
are not included in the fit, because they lead to significantly increased ;^^/dof in the fit and hence 
we conclude that NLO ;^PT is not appropriate for such large mass-values. 



Our data for /ps and mps is described excellently by the ;^PT formulae, see figures 5(b), 6(a) 



and 6(b), where we plot appropriate, dimensionless ratios. The fitted values of the parameters are 

A3//o = 6.41(26), A4//o = 11.51(21), 
«/o 1/3=3.9 = 0.0527(4) , a/o 1/3=4.05 = 0.0411(4) , (4.9) 
2aZ^Bo|/3=3.9 = 4.87(4) , 2aZ^Bo|/3 =4.05 = 3.76(3) , 

with ;f^/dof =12/12 and statistical errors only. The sensitivity to A3_4 is visualised by the deviation 



from linearity in figures 5(b) and 6(b). 
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The statistical uncertainties are estimated using a bootstrap procedure: at each of our data 
points we produced 1000 bootstrap samples of /ps and mps and used them to perform 1000 fits. 
The statistical uncertainty is then given by the variance over these 1000 fit results. Note that this 
procedure automatically takes into account the cross-correlation between /ps and mps for a given 
ensemble. 

Systematic uncertainties are - as usual - harder to estimate. We include the effects coming 
from (a) finite size effects by using GL instead of CDH to estimate finite volume effects, (b) finite 
size effects by including A1//0 and A2//0 as free parameters in the fit and (c) lattice artifacts by 
performing the fits separately for j3 = 3.9 and j8 = 4.05. They are discussed in more detail in 



Ref. [44] 



Lattice Calibration and Low Energy Constants 

As mentioned above, by fixing the value of /ps to the physical value of the pion decay constant 
fn = 130.7 MeV where nips/fps assumes its physical value we can calibrate our lattices for the two 
/3 -values, with the result 

a|^^3 9 =0.0855(5) fm, a||3=4.05 = 0.0667(5) fm. (4.10) 

The ratio a|/j=3.9/a||3=4.o5 = 1-28(1) can be compared with the ratio ('"o/«)|j3=4.o5/(''o/^i')|/5=3.9 = 
1.27(1) and we find excellent agreement. Moreover, using those values for the lattice spacings we 
can get an estimate of the ratio Z^|jg=3.9/-Z^ 1/3=4.05 from our fit and compare it to the correspond- 
ing ratio as determined with RI-MOM. Also here the agreement is excellent. We take this as an 
indication that the combined fit does not hide lattice artifacts in fit parameters. 
Having set the scale allows us to determine the low energy constants 

£3 = log(4) =3.44(8)(35), £4 = log =4.61 (4) (11), (4.11) 

where we set nij^ = 139.6 MeV conventionally. The first error is statistical, the second systematical 



(see Ref. [ 144| ] for details on how we estimate the systematic errors). The statistical accuracy on 



these estimates for £3 4 is quite impressive and the mean values are in good agreement with the 



literature (see Refs. [51, 52] for recent reviews). 



We conclude by mentioning that the values of £3 4 determined from the aforementioned ;^PT fit 



for /ps and mps in the continuum [ |44| ] are in good agreement with the values quoted in Eq. ( [4. lit ). 
4.2 NucleonMass 

In figure [7(a)| we plot our preliminary results for the nucleon mass in units of ro as a 
function of (romps )^. Also in this quantity scaling violations are not visible within the statistical 
accuracy. Moreover, as shown by the comparison of ensemble fii and Bd, finite size effects seem 
to be neghgible. Concentrating on jS = 3.9 only, since we do not yet have enough data points at 



j3 = 4.05, we fit the leading one-loop Heavy Baryon ;^PT formula [53, 54] 



= mg, + 4cim^s " '"ps (4.12) 

3271/^5 
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Figure 7: (a) Nucleon mass rgmN as a function of (romps)^ for j3 = 3.9 and j3 = 4.05. (b) (''o(('«ps)^ 
(mpg)^) as function of (a/ro)^ for ensembles B\ and C\. 



with gA = 1.267 and the two free parameters and c\ to our data. The fit provides a good 
description of the data with ;^^/dof = 0.2. With the physical nucleon mass as input the value of the 
lattice spacing at j8 = 3.9 can be determined to a = 0.0879(12) fm in very good agreement with 
the determination from fj^. This result not only successfully cross checks the determination of the 
lattice spacing in physical units from fj^. If confirmed, it also represents a test of QCD as the theory 
of strong interactions. For more details, including the error determination and other baryon masses 



see Ref. [55|. 



4.3 Flavour Symmetry Breaking Effects 

As mentioned in section ^ flavour symmetry is explicitly broken. As a consequence there is 
a potential difference of ff{a^) between the masses of charged and neutral pseudo scalar mesons. 
Notice that to the latter also disconnected diagrams contribute. We have determined the mass of 
the neutral pseudo scalar meson for the ensembles Bi, B3, B^, C\ and C2. The results are shown in 
table ^ where we also report the the corresponding values for the charged meson mass. 

The neutral pseudo scalar meson is lighter than the charged one. This observation is consis- 



tent with ;^PT predictions for the observed first order phase transition scenario (see Ref. [ |13| ] and 
references therein). For ensembles Bi and C\ we plot in figure 7(b)| the difference 



4{H^f-{m^,f) (4.13) 



as a function of (ro/a) . The quantity (4.13) is expected to scale linearly in {ro/a) towards the 



continuum, which is confirmed by our data. The dashed line in figure 7(b)| is not a fit, but it is there 
only to guide the eye. 

Even if this analysis provides evidence for the expected scaling behaviour of the pion mass 
splitting, the effect still amounts to about 16% at our smallest value of the lattice spacing and 
charged wips ~ 300 MeV. If compared to the "natural" size one would expect for ff{a^) cut-off 
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Ensemble 


arrips 


anips 


amy 


amy 


Bi 


0.13623(65) 


0.111(6) 


0.404(22) 


0.395(17) 


B3 


0.19403(50) 


0.167(9) 


0.428(08) 


0.419(17) 


Be 


0.13377(24) 


0.110(8) 


0.416(12) 


0.400(25) 


Ci 
Ci 


0.1038(6) 
0.1432(6) 


0.091(11) 
0.126(06) 


0.337(20) 
0.337(12) 


0.372(29) 
0.346(12) 



Table 3: Comparison of values for the charged and the neutral pseudo scalar (vector) meson masses. 



effects, which is of the order of a^Agf^^, one has to conclude that this effect is unexpectedly large 
and one would like to be able to better understand it. Of particular interest is the question whether 
this effect arises dynamically, or is due to large coefficients in the Symanzik expansion. In the latter 
case one might expect that this effect is not restricted to only the neutral pseudo scalar meson mass. 

Note that all other possible splittings we have determined so far are negligible. For instance 
the splittings in the vector and A channels appear to be consistent with zero (see table ^ for my and 



Ref. 1 55] for the A). Similarly, the difference between the decay constants of charged and neutral 
pseudo scalar meson is negligible. 

Currently, the ETM collaboration is investigating this question [^, 57 1, and from a theoretical 



point of view, an analysis a la Symanzik of the charged and the neutral pseudo scalar meson masses 
lead to the formulae 



{mlsf = ml + a^^j: + £ 
{m^sf =ml + ff{a^ml. 



[a^m^ja"^) 



(4.14) 



which show that the difference (nip^)^ — (nip^)^ is given by the term proportional to l^ji. Here 
Cji = (^"l-S^el^r") and is the dimension six term in the Symanzik effective Lagrangian. 

The main result of the analysis of Refs. |5^ is that Q is a large number which in the vac- 
uum saturation approximation can be estimated to be proportional to \Gk\^, where Gji = {0\P^\k'^). 
The latter matrix element is numerically large: one finds |G;r|7AQCD around 20- 25 [|6[|7|]. This 
result can provide an interesting physical explanation for the large splitting observed in the pseudo 
scalar masses. Moreover, since it can be shown that l^ji enters only the neutral pseudo scalar mass 
(and related quantities), one also finds a possible explanation of why all other splittings determined 
so far turn out to be small. In addition it provides hope that in the future we shall not find large 
lattice artifacts due to flavour breaking. 

At this point the question might arise whether the large splitting in the pseudo scalar masses 
affects the ;^PT fits and in particular the ^PT FS estimates. At the current level of our analysis, 
where we work under the assumption that lattice artifacts are zero in charged mps and /ps - and 
we have shown good theoretical and numerical evidence for this - all the fits are performed in the 
continuum, where flavour symmetry is restored. Hence, we do not expect any effect of flavour 
symmetry breaking on our analysis with ;tPT. 

4.4 Further ETMC Results 

Our current simulations contain two light quark flavours with degenerate mass in the sea. In 
the unitary set-up we can hence determine the average up-down quark mass in the MS scheme by 
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using x^T fits discussed in section pTl] The result is m^^{jj.R = 2 GeV) = 3.62(10)(23) MeV. We 
have again used the estimate for Zp coming from Rl-MOM |2^]. Here and in the following the 
first error is statistical, the second systematical. 

In order to also determine the strange and possibly the charm quark masses, we have to use 
a partially quenched set-up where we compute propagators on the available gauge configurations 
with several values of the valence quark mass - which are now different from the sea quark mass 
- around the strange and the charm quark mass. The set-up and the calculation are described in 
detail in Refs. ^ and the result for the strange quark mass for j8 = 3.9 only reads m^^{iJ.R = 
2 GeV) = 105(3) (8) MeV. A comparison to other lattice QCD determinations of mg shows that 
non-perturbative renormalisation has significant impact on the final result. Estimates obtained with 
perturbative renormalisation tend to be significantly lower than estimates with non-perturbatively 
computed renormalisation constants. For more details as well as for further related results, such as 
the kaon decay constant /k and the ratio /k/Ztt. we refer to Ref. [59, ^ ]. 

In Ref. [|60|] we present our current estimates for the charm quark mass at jS = 3.9. Our result 



reads m^^ (^Ur = m^^) = 1 .481 (21) (94) GeV. In Ref. [|o|| also results for fo and fo, are presented, 
which are in good agreement with experiment. 

In addition, the gauge configurations at jS =3.9 have been used to calculate a variety of 3-point 
correlation functions relevant for semileptonic weak decays and electromagnetic transitions of light 
and heavy-light pseudo scalar mesons. Results for the the vector, scalar and tensor form factors of 
the pion, the vector and scalar form factors relevant for K13 decay and for the Isgur-Wise function 
in the heavy-quark limit can be found in Ref. [51 1. On the same set of gauge configurations there 
are preliminary results for the first moment of the pion quark distribution function available. They 
can be found in [62] together with an effective stochastic method to determine this quantity in 
lattice simulations. 

The ETM collaboration has made a serious attempt to determine properties of flavour singlet 
mesons with Mtm-LQCD, a first account of which can be found in Ref. [63]. The most interesting 
result of this investigation is that the mass of the 772 meson (the one related to the anomaly in two 
flavour QCD, i.e. not a Goldstone boson) is consistent with a constant behaviour in the chiral limit 
and the value around 800 MeV is compatible with expectations from a model computation []64|]. 
(See also Ref. [||]) 

For certain quantities, like Bk, it is crucial that the lattice formulation exhibits good chiral 
properties. The overlap operator, which obeys exact chiral symmetry at finite value of the lattice 
spacing, can be used in a mixed action approach as a valence operator on a maximally twisted mass 
sea. An exploratory study and first results profiting from the chiral properties in the valence sector 
can be found in Ref. [|6^. 

A systematic uncertainty we cannot control at the moment is due to the fact that the effects of 
the strange quark are not taken into account in the simulations. In order to include these effects in 
Mtm-LQCD, maintaining ff{a) improvement at the same time, a 1 + 1 split heavy doublet of quarks 
has to be simulated in addition to the mass degenerate light quark doublet [|67||. In an exploratory 



study, published in Ref. [58], it was shown that this approach is feasible and that tuning is possible. 
However, it was also shown that the effect of the first order phase transition, mentioned in previous 
sections, strengthens significantly when the heavy doublet is added. In Ref. [^ we report on an 
attempt to cure this potential problem by using stout smearing [[70[]. Our preliminary results suggest 
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that stout smearing reduces indeed the effects of the phase transition significantly. 

An interesting investigation at tree level of perturbation theory is presented in Ref. [f/TI]. In 
this framework cut-off effects can be studied using analytic calculations. Scaling properties of 
Wilson and Wilson twisted mass fermions are compared and it is shown that maximally twisted 
mass fermions scale with a rate of 



For a study using Mtm-LQCD for simulations of QCD thermodynamics see Ref. []72||. 

4.5 Exploring the £ -regime with Maximally Twisted Mass Fermions 

A slightly different direction as compared to all the aforementioned results is explored by 
performing studies of the £-regime [ [7^ ] with Mtm-LQCD at j8 = 3.9 [74]. Simulations in the £- 



regime are not restricted to formulations with exact chiral symmetry. And, since the twisted mass 
parameter provides - as discussed before - an infra-red cut-off to the eigenvalue spectrum of the 
twisted mass Dirac operator there is also no technical complication to expect, unlike that arises e.g. 
for Wilson fermions when the quark mass becomes too small. 



Preliminary results of this investigation presented in Ref. [ 74 ] are quite encouraging. Simula- 
tions turn out to be feasible, perform smoothly and are much less computer time demanding than 
simulations with e.g. the overlap operator. A first result of this study at j3 = 3.9 is an estimate for 
the chiral condensate ^{\Ir = 2 GeV) = (264(12) (4)lgo MeV)^ which is in perfect agreement 
with the value determined from the ^PT fit 1^{ixr = 1 GeV) = (267(4) (7) MeV)\ 

5. Summary 

In this proceeding contribution we have summarised the current status of the two flavour 
project of the European Twisted Mass collaboration. The collaboration has generated gauge con- 
figurations using a doublet of mass degenerate Wilson twisted mass fermions at maximal twist for 
three different values of the lattice spacing a < 0.1 fm, volumes with physical extent larger than 
2 fm and values for the (charged) pseudo scalar meson mass in the range of 300 to 700 MeV. The 
results for the two finer lattice spacings of about a ~ 0.086 fm and a ~ 0.067 fm are close to final, 
while the results for the coarsest lattice spacing are still in a preliminary state. 

We have presented results for (charged) /ps and mps and a scaling analysis for /ps at fixed, but 
finite volume indicating that ff{a) improvement in Mtm-LQCD works very well. Lattice artifacts 
turn out to be compatible with zero to our current statistical accuracy, which is of the order of 1%. 

We have shown evidence that finite size effects in /ps and mps can be described by means of 
formulae derived in chiral perturbation theory to the level of statistical accuracy of the data. It turns 



out that the asymptotic Luscher formula presented in Ref. Q49I ] works better than the NLO formula 



from Ref. [46] 



NLO continuum chiral perturbation theory [ |46| ] can be used successfully to describe the quark 



mass dependence of mass and decay constant of the charged pseudo scalar meson. The low energy 
constants Bq, Fq, It, and £^ can be determined to high statistical accuracy. The corresponding fit can 
also be used to determine the lattice spacings using the physical values of fj^ and rrij^. The results 
are in good agreement with the determination from the nucleon mass. 

There is theoretical and numerical evidence that large flavour breaking effects appear only in 
the mass of the neutral pseudo scalar meson and trivially related quantities. In particular, all other 
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flavour splittings measured so far turn out to be compatible with zero. In addition the mass splitting 
in the pseudo scalar meson masses scales as expected towards the continuum limit. 

The collaboration is analysing the available gauge configurations for many more physical 
quantities. As an example we have presented first results for the strange and the charm quark 
masses. For the strange quark mass, where we can compare to other lattice determinations, it turns 
out that the difference between perturbative and non-perturbative renormalisation is significant. 
Quantities like /k/A or /d//d, are in good agreement with experiment. 

In the future we plan to repeat all these calculations with strange and charm quark effects 
taken into account in the simulations. The ETM collaboration is currently investigating the optimal 
set-up for simulations with 2+1 + 1 dynamical quark flavours and maximal twist. Algorithms and 
codes are available and the simulations are due to start. 

We conclude by mentioning that all ETMC ensembles are stored on ILDG disk space [[7^. 
They are available to non-ETMC members on a request basis, whenever there is no overlap to 
ongoing ETMC projects. 
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